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(Dated: March 24, 2015) 

Using data samples collected with the BESIII detector operating at the BEPCII collider at center- 
of-mass energies from 3.810 to 4.600 GeV, we perform a study of e''"e“ vJ/'^ s-nd TT^Jlif. Sta¬ 

tistically significant signals of e^e“ —>■ 'gj/ip are observed at ^ys = 4.190, 4.210, 4.220, 4.230, 
4.245, 4.260, 4.360 and 4.420 GeV, while no signals of e^e“ —>■ Tv^J/ip are observed. The measured 
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energy-dependent Born cross section for e'*"e“ —>■ rjj/tp shows an enhancement around 4.2 GeV. The 
measurement is compatible with an earlier measurement by Belle, but with a significantly improved 
precision. 

PACS numbers: 13.25.Gv, 13.66.Bc, 14.40.Pq, 14.40.Rt 


I. INTRODUCTION 

During the last decade, new charmoniumlike vector 
states, such as the F(4260), F(4360) and y(4660), 
have been observed by BABAR [1, 2], Belle [3-6] and 
CLEO [7]. The masses of these new Y states are above 
the DD production threshold, ranging from 4.0 to 4.7 
GeV/c^. Since all of them are produced in e“'"e“ anni¬ 
hilation (either directly or via the initial state radiation 
(ISR) process), and since they have been observed to de¬ 
cay in dipion hadronic transitions to the J/tp or ■0(3686), 
one would naturally interpret these states as vector char- 
monium excitations. However, peculiar features of these 
Y states reveal an exotic nature that likely excludes a 
conventional charmonium interpretation. These features 
include a discrepancy with the spectrum of vector char¬ 
monium states predicted by the potential model given in 
reference [8], a surprisingly large coupling to final states 
without open-charm mesons [9, 10], and a lack of obser¬ 
vation in the inclusive hadronic cross section [11]. Also, 
very recently, several charged charmoniumlike states — 
the Zc(3900)± [4, 12, 13], Zc(3885)± [14], Zc(4020)± [15], 
Zc(4025)^ [16], as well as their isospin partners, the neu¬ 
tral states Zc(3900)° [13] and Zc(4020)° [17] — were ob¬ 
served in the same mass region as these Y states. This 
suggests that the nature of the Y states could be re¬ 
lated to that of the Zc states. Moreover, BESIII re¬ 
cently reported on the measurement of the cross sec¬ 
tions of e+e” —>■ 7r+7r“/ic [15] and e+e" wyco [18]. 
The observed cross sections as a function of center-of- 
mass (CM) energy are inconsistent with the line shape 
of e+e“ —>■ 7r+7r“ J/0 [1]. These observations hint at the 
existence of a more complicated and mysterious underly¬ 
ing dynamics. 

Many theoretical interpretations have been proposed 
to classify these Y states, such as hybrid charmo¬ 
nium [19], tetraquark [20], or hadronic molecule [21] mod¬ 
els, but none of them has been able to describe all ex¬ 
perimental observations in all aspects. Searching for new 
decay modes and measuring the line shapes of their pro¬ 
duction cross sections will be very helpful for these Y 
states interpretation. Hadronic transitions (by rj, 7r°, or 
a pion pair) to lower charmonia like the J/ip are also 
regarded as sensitive probes to study the properties of 
these Y states [22]. 

The cross sections of e+e“ —>■ fjj/'ip and tt°J ftp above 
the DD production threshold have been evaluated within 
a non-relativistic framework derived from QCD [23] , and 
their line shapes are predicted to be strongly affected by 
open charm effects [24]. Belle, BESIII and CLEO-c have 
measured the production cross sections of e+e“ —?► ??T/0 
above the open charm threshold [25-27]. However, Belle 


and CLEO-c results suffer from large statistical uncer¬ 
tainties. BESIII reported on a more accurate result, but 
the measurement was limited to a single center of mass 
energy of ^/s = 4.009 GeV. Experimental studies with 
large data samples in a broad energy region may shed 
light on the nature of the Y states. 

In this paper, we report a measurement of the Born 
cross sections of e“'"e“ —>■ riJ/ip and ■K^J/'p from e“'"e“ 
CM energies 3.810 GeV to 4.600 GeV with data samples 
taken by BESIII. In our analysis, the rj and 7r° are recon¬ 
structed in their two-photon decay mode and the J/ip via 
its decay into lepton pairs (£+£“). 


II. DETECTOR AND MONTE CARLO 
SIMULATION 

BEPCH [28] is a double-ring e+e” collider running at 
CM energies ranging from 2.0 to 4.6 GeV, and providing a 
peak luminosity of 0.85 x 10^^ cm“^s“^ at the CM energy 
of 3.770 GeV. The BESIII [28] detector has a geometrical 
acceptance of 93% of 47r and has four main components. 
(1) A small-cell, helium-based (40% He, 60% CaHg) main 
drift chamber (MDC) with 43 layers provides an average 
single-hit resolution of 135 ^m, and a charged-particle 
momentum resolution in a 1 T magnetic field of 0.5% at 
1.0 GeV/c. (2) A time-of-flight system (TOE) is con¬ 
structed of 5 cm thick plastic scintillators, with 176 de¬ 
tectors of 2.4 m length in two layers in the barrel and 
96 fan-shaped detectors in the endcaps. The barrel (end- 
cap) time resolution of 80 ps (110 ps) provides a 2cr Kj-K 
separation for momenta up to GeVjc. (3) An elec¬ 

tromagnetic calorimeter (EMC) consists of 6240 CsI(Tl) 
crystals in a cylindrical structure (barrel) and two end- 
caps. The photon energy resolution at 1.0 GeV/c is 2.5% 
(5%) in the barrel (endcaps), and the position resolution 
is 6 mm (9 mm) in the barrel (endcaps). (4) The muon 
system (MUG) is located in the iron flux return yoke 
of the superconducting solenoid and consists of 1000 m^ 
of Resistive Plate Chambers (RPCs) in nine barrel and 
eight endcap layers. It provides a position resolution of 
2.0 cm. 

The optimization of the selection criteria, the deter¬ 
mination of detection efficiencies, and the estimations 
of potential backgrounds are performed based on Monte 
Carlo (MC) simulations taking various aspects of the ex¬ 
perimental setup into account. GEANT4-based [29] MC 
simulation software, which includes the geometric and 
material description of the BESIII detector, the detec¬ 
tor response and digitization models, as well as an ac¬ 
counting of the detector running conditions and perfor¬ 
mances, is used to generate MC samples. In the Simula- 
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tion, the electron-positron collision is simulated with the 
KKMC [30, 31] generator taking into consideration the 
spread in the beam energy and ISR. In this analysis, large 
signal MC samples of e+e“ —?► s-nd e+e“ —?► TT^J/ip 

are generated at CM energies corresponding to the ex¬ 
perimental values, where the line shape of the produc¬ 
tion cross section of these two processes, assumed to be 
identical, are taken from the Belle experiment [25]. 


III. EVENT SELECTION 

The candidate events of e+e” —>■ riJ/tp and 
are required to have two charged tracks with a total net 
charge of zero and at least two photon candidates. 

Charged tracks are reconstructed from the hits in the 
MDC. Each charged track is required to have a polar an¬ 
gle that is well within the fiducial volume of the MDC, 

I cos0| < 0.93, where 0 is the polar angle of the track in 
the laboratory frame, to have a point of closest approach 
to the interaction point that is within ±10 cm along the 
beam direction and within 1 cm in the radial direction, 
and to have a momentum p larger than 1.0 GeV/c. Elec¬ 
tron and muon separation is carried out by making use 
of the deposited energy in the EMC. Tracks with an en¬ 
ergy deposition oi E < 0.4 GeV are identified as muons, 
while tracks with E/p > 0.8 c are identified as electrons 
or positrons. 

Photon candidates are reconstructed by isolated show¬ 
ers in the EMC. The photon energy is required to be at 
least 25 MeV in the barrel (| cos 9\ < 0.80) and 50 MeV in 
the endcaps (0.86 < |cos0| < 0.92). To eliminate show¬ 
ers produced by charged particles, the angle between the 
shower and the nearest charged track must be larger than 
20 degrees. To suppress electronic noise and energy de¬ 
positions unrelated to the physical event, the EMC time 
t of the photon candidate must be in coincidence with 
collision events, in the range from 0 < t < 700 ns. 

A kinematic fit that imposes momentum and energy 
conservation (4C) is implemented under the hypothesis 
of e“'"e“ —> to improve the momentum and en¬ 

ergy resolutions of the final-state particles and to reduce 
the potential backgrounds. The chi-square of the kine¬ 
matic fit, is required to be less than 40. If there 
are more than two photons in an event, the combina¬ 
tion of with the least is chosen. To sup¬ 

press the backgrounds from radiative Bhabha and ra¬ 
diative dimuon events associated with a random photon 
candidate, the energy of each selected photon is further 
required to be larger than 80 MeV. 

Figure 1 depicts scatter plots of the invariant mass 
of lepton pairs, M(£+£“), versus that of two photons, 
M( 77 ), for data taken at %/s = 4.230 and 4.260 GeV. A 
clear accumulation of events is observed around the in¬ 
tersection of the r\ and J/i/ mass regions, which indicates 
e+e” r]J/ -Ip signals. There is no significant signal ob¬ 
served around the intersection of the and J/'ip mass 
regions. MC studies show that dominant backgrounds 


are from the radiative Bhabha and dimuon events, and 
are expected to be distributed uniformly around the J/ip 
and p/t:^ mass regions. A significantly larger background 
yield is observed in the e“'"e“ mode than in the 
mode, which is due to the much larger Bhabha scatter¬ 
ing cross section compared with the dimuon cross sec¬ 
tion. MC simulations show that the resolution of the 
invariant mass distributions of lepton pairs is about 10.7 
MeV/c^ for the pT mode and 11.5 MeV/c^ for the 
e+e" mode. The candidate event of e+e” —>■ pJ/ip 
is required to be within the J/ip signal region, defined 
as 3.067 < < 3.127 GeV/c^. Sideband re¬ 

gions, defined as 2.932 < M{£~^i~) < 3.052 GeV/c^ and 
3.142 < < 3.262 GeV/c^, four times as wide 

as the signal region, are used to estimate the non-J/'0 
background contributions. 

After selecting the J/ip signal, the invariant mass dis¬ 
tributions of two photons, M( 77 ), are shown in Fig. 2 for 
data at ^/s = 4.230 and 4.260 GeV. Clear p signals are 
observed. The corresponding normalized distributions 
from the events in the J/ip sideband regions are shown 
as shaded histograms in the plots. The backgrounds are 
well described hy J/ip sideband events and show no peak¬ 
ing structure within the p signal region. 

The process e+e“ —>■ ir^J/ip is also searched for in 
the J/ip —> tJ-~ mode by analyzing the dis¬ 

tribution around the mass region. Such a search is 
not performed for the J/ip e^e~ mode due to the 
large background of radiative Bhabha events. Due to the 
misidentification of as peaking background from 
e+e” —>■ 7 r+ 7 r“ 7 r° would contaminate the 7 r° signal for 
both candidate events within the J/ip signal or sideband 
regions. To remove such backgrounds, we require that 
at least one charged track has a muon counter hit depth 
larger than 30 cm. Figure 3 shows the M{yy) distribu¬ 
tions around the mass region after this requirement. 
No significant signal is observed for 7 r° —>■ 77 decays. 


IV. FITS TO THE M(77) SPECTRUM AND 
CROSS SECTION RESULTS 

After imposing the J/ip signal selection, an unbinned 
maximum likelihood fit is performed on ± 7 ( 77 ) in the 
J/ip ^ e+e“ and pi~ modes, respectively. The prob¬ 
ability density function (PDF) of the M{yy) distribu¬ 
tion for p signals is obtained from signal MC simulations 
convoluted with a Gaussian function, where the Gaus¬ 
sian function describes the difference in resolution be¬ 
tween data and MC simulation, and its parameters are 
left free in the fit. The background shape is described by 
a second-order Chebyshev polynomial function. The cor¬ 
responding fit results for ^/s = 4.230 and 4.260 GeV are 
shown in Fig. 2 and the numbers of p signal events are 
summarized in Table I. The statistical significances of p 
signals are larger than 8 (t, which are examined using the 
differences in log-likelihood values of fits with or without 
an p signal component. 
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FIG. 1. Scatter plots of ) versus Mi^'yy) for data at yfs = 4.230 (top panels (a, b)) and 4.260 GeV (bottom panels (c, 

d)). The two panels on the left-hand side correspond to the mode and the right-hand side to the e''"e“ mode. The blue 

dotted lines denote the rf/'K^ and J/'i/j mass bands. The red dashed lines denote the sideband regions of J/iji. 


TABLE I. Results on e'*'e“ in data samples in which a signal is observed with a statistical significance larger than 

5(t. The table shows the CM energy y/s, integrated luminosity £int, number of observed rj events 

from the fit, efficiency radiative correction factor (1 -I- 5’^), vacuum polarization factor (1 -|- 5”), Born cross section 

fj,~)/a^{e^e~) and combined Born cross section crcom- Ths first uncertainties are statistical and the second systematic. 


Ui(GeV) jC(pb ^ 

) N°°-{e+e-) 

Em(%) 

ee(%) (1±(SG (1±5U 

( 71+77 )(pb) 

cr"(e+e )(pb) 

<^Com(pfi) 

4.190 

43.1 

17.5 ±4.3 

10.4 ±3.6 

35.2 

24.1 

0.866 

1.056 

53.7 ± 13.2 ±3.1 

46.6 ±16.1 ±1.7 

50.8 ±10.2 ±2.1 

4.210 

54.6 

25.7 ±5.1 

14.8 ±4.5 

33.7 

23.1 

0.914 

1.057 

61.6 ± 12.2 ±4.1 

51.7 ± 15.7 ±4.5 

57.8 ±9.6 ±3.2 

4.220 

54.1 

32.6 ±5.8 

11.4 ±3.9 

33.1 

22.8 

0.937 

1.057 

78.2 ± 13.9 ±5.0 

39.6 ± 13.6 ±2.9 

57.7 ±9.7 ±3.0 

4.230 

1091.7 

394.3 ± 20.9 

274.9 ±20.1 

32.4 

22.3 

0.960 

1.056 

46.8 ±2.5 ±2.5 

47.3 ±3.5 ±3.4 

47.0 ±2.0 ±2.2 

4.245 

55.6 

9.3 ±3.3 

9.7 ±3.6 

31.4 

21.7 

0.992 

1.056 

21.6 ±7.7 ±2.4 

32.6 ±12.1 ±3.5 

24.8 ±6.5 ±2.0 

4.260 

825.7 

94.4 ± 10.5 

75.9 ± 11.9 

30.3 

20.9 

1.021 

1.054 

14.9 ±1.7 ±1.1 

17.4 ± 2.7 ± 1.2 

15.7 ± 1.4 ±0.9 

4.360 

539.8 

19.8 ±5.3 

23.9 ±7.7 

25.7 

17.7 

1.168 

1.051 

4.9 ± 1.3 ±0.7 

8.7 ± 2.8 ± 1.1 

5.6 ±1.2 ±0.6 

4.420 

1074.7 

56.9 ±8.2 

42.6 ±9.9 

24.2 

16.7 

1.225 

1.053 

7.3 ±1.1 ±0.8 

7.8 ±1.8 ±0.7 

7.5 ±0.9 ±0.6 


The same event selection criteria are implemented on 
the other 15 data samples taken at different CM energies. 
We observe a significant yield of ry signal with more than 
5(7 statistical significance for data at y/s = 4.190, 4.210, 
4.220, 4.245, 4.360 and 4.420 GeV. Using the same fit pro¬ 
cedure, the numbers of 77 signal events for these energies 
are also obtained and listed in Table I. There are no sig¬ 
nificant 77 signals observed for the other 9 energy points 
(Vs = 3.810, 3.900, 4.090, 4.310, 4.390, 4.470, 4.530, 
4.575, 4.600 GeV), and the upper limits at the 90% con¬ 
fidence level (G.L.) on the Born cross section are deter¬ 
mined with the J/tp ^ decay mode only, due to the 

large background from Bhabha events in J/tp —^ e+e“. 


Since the statistics are low for the no-signal data samples, 
the number of observed events is obtained by counting 
the entries in the 77 signal region (0.518 < < 0.578 

GeV/c^). The number of background events in the signal 
region is estimated by the events in the 77 sideband region 
or J/ip sideband region (with an additional 77 signal mass 
window requirement) by assuming a flat distribution of 
background around signal regions. The 77 sideband re¬ 
gion is defined as 0.383 < < 0.503 GeV/c^ and 

0.593 < M{jj) < 0.713 GeV/c^, where their sizes are 
four times as that of the signal region. The results are 
all listed in Table II. 
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FIG. 2. Invariant mass distributions of two photons for data at y/s = 4.230 (top panels (a, b)) and 4.260 GeV (bottom panels 
(c, d)). The left two plots are for the J/^|) fi~ mode and the right two for J/?/? e'^e~ mode. Dots with error bars are 

for data in the J/^p signal region, the green shaded histograms for the normalized J/xjj sideband events, the red solid curves for 
the total fit results and the blue dotted curves for the background from the fit. 
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FIG. 3. Invariant mass distributions of two photons for data at y/s = 4.230 (a) and 4.260 GeV (b) in J/'tp —>■ fj,^ mode. Dots 
with error bars are for data in the J/ij) signal region, the green shaded histograms are the normalized J/'i/j sideband events and 
the red histograms are 7r° MG signal with arbitrary normalization. 


The Born cross section is calculated by: 

Anf + + 

where N°'°^ is the number of observed signal events, £int 
is the integrated luminosity, (I+ (5’') is the ISR correction 
factor which is obtained by QED calculation [32] and 
taking the line shape of the Born cross section measured 


by the Belle experiment [25]. The vacuum polarization 
(VP) factor (1 + 5'') is taken from a QED calculation 
with an accuracy of 0.5% [33], e is the detection efficiency 
including reconstruction and all selection criteria, B is the 
product branching ratio, and B{J/'ip —>■ £+£“) • B{ri —>• 
77 ), taken from the Particle Data Group (PDG) [34]. 

The hnal Born cross sections of e+e” —>■ rjj/i!) at en¬ 
ergy points with a statistically significant observation of 
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TABLE II. Upper limits of e^e“ 'nJ/'4’ using the fJ-~ mode. The table shows the CM energy ^/s, integrated luminosity 
Tint, number of observed r] events number of background from 77 sideband N^, and from J/ip sideband efficiency e, 

upper limit of signal number with the consideration of selection efficiency /e. (at the 90% C.L.), radiative correction factor 
(1 + (5’’), vacuum polarization factor (1 + <5”), Born cross section ct® and upper limit on the Born cross sections a^p (at the 90% 
C.L.). The first uncertainties are statistical and the second systematic. 


yi(GeV) 

C{ph-^) 

crq 

jysb 

lyjj 


e(%) 


(1 + V) 

(1-f <5”) 

a“(pb) 

o-up(pb) 

3.810 

50.5 

5 

9 

11 

32.3 

< 

23.3 

1.243 

1.056 

5.501:^ ±0.2 

< 15.1 

3.900 

52.6 

5 

8 

7 

38.3 

< 

20.9 

0.775 

1.049 

7.9l^ g ±0.3 

< 20.8 

4.090 

52.6 

7 

7 

5 

31.0 

< 

36.3 

1.087 

1.052 

12.27-§ ±0.9 

< 25.9 

4.310 

44.9 

1 

4 

2 

27.4 

< 

11.7 

1.105 

1.053 

0.0029 ±0.0 

< 9.5 

4.390 

55.2 

5 

1 

4 

25.0 

< 

38.4 

1.198 

1.051 

11.705 4 ±0.6 

< 23.5 

4.470 

109.9 

2 

12 

8 

23.5 

< 

14.7 

1.258 

1.055 

-1.200® ±0.1 

< 4.3 

4.530 

110.0 

5 

6 

4 

22.8 

< 

38.2 

1.295 

1.055 

4.3001 ± 0.2 

< 10.8 

4.575 

47.7 

2 

2 

1 

22.6 

< 

22.5 

1.314 

1.055 

4.300^ ± 0.2 

< 14.5 

4.600 

570.0 

5 

34 

19 

22.4 

< 

23.0 

1.323 

1.055 

-0.800^ ±0.1 

< 1.2 


signal events are listed in Table I. 

For the other energy points where the rj signal is not 
significant, we set upper limits at the 90% C.L. on the 
Born cross section. The upper limit is calculated by a 
frequentist method with a profile likelihood treatment 
of systematic uncertainties taken into account in the effi¬ 
ciency uncertainty, which is implemented by a C-|—h class 
TROLKE in the ROOT framework [35]. The numbers of 
observed signal events and estimated background events 
are assumed to follow a Poisson distribution, and the effi¬ 
ciencies are assumed have Gaussian uncertainties. Since 
the number of background events can be estimated from 
either the rj or J/ip sideband events, the one with the 
larger upper limit on the Born cross section is taken as 
the final result as a conservative estimation. The results 
on the upper limits are listed in Table II. 

Since there is no significant signal of e+e” —>■ 
observed at any energy, we set upper limits at the 90% 
C.L. on the Born cross section. The number of observed 
events is obtained by counting the entries in the tt*^ signal 
region (0.120 < M( 77 ) < 0.150 GeV/c^). The number 
of background events in the signal region is estimated 
by counting the number of events in the tt^ sideband 
regions (0.055 < M{jj) < 0.115 GeV/c^ and 0.155 < 
M( 77 ) < 0.215 GeV/c^) or J/t/i sideband regions (with 
an additional signal mass window requirement). The 
same frequentist method is implemented to extract the 
upper limits. The results and the related variables used 
to calculate the upper limit are listed in Table III. 


V. SYSTEMATIC UNCERTAINTIES 

Several sources of systematic uncertainties are con¬ 
sidered in the measurement of the Born cross sections. 
These include differences between data and MC simula¬ 
tion for the tracking efficiency, photon detection, kine¬ 
matic fit, mass window requirement, the fit procedure, 
the shower depth in the MUG, MC simulation of the ISR 
correction factor and vacuum polarization factor, as well 


as uncertainties in the branching fractions of intermedi¬ 
ate state decays and in the luminosity measurements. 

(a) Tracking: The uncertainty of the tracking effi¬ 
ciency is investigated using a control sample '0(3686) —>■ 

with the subsequent decay of J /'0 —>■ . 

The difference in tracking efficiency for the lepton recon¬ 
struction between data and MC simulation is estimated 
to be 1% per track. So, 2% is taken as the systematic 
uncertainty for the two leptons. 

(b) Photon detection efficiency: The uncertainty due 
to the photon detection and reconstruction efficiency is 
1% per photon [36]. This value is determined from stud¬ 
ies using background-free control samples J/tp —^ / 9 ° 7 r° 
and e+e” —>■ 77 . Therefore, an uncertainty of 2% is 
taken for the detection efficiency of two photons. 

(c) Kinematic fit: In order to reduce the difference on 
the 4C kinematic fit xlc between data and MC simula¬ 
tions, the track helix parameters (0o, tan A) of sim¬ 
ulated tracks have been corrected, where 00 is th® az¬ 
imuthal angle that specifies the pivot with respect to the 
helix center, k is the reciprocal of the transverse momen¬ 
tum, and tan A is the slope of the track. The correction 
factors are obtained from a nearly background-free sam¬ 
ple of e+e” —>■ -TT+Tr” J/'0 and J/ip —>■ e'^e~ jiffiyT at 

= 4.230 GeV. An alternative detector efficiency is 
evaluated with the same MC samples, but without helix 
parameters corrections. The difference in this efficiency 
from its nominal value is taken to be the uncertainty due 
to the 4C kinematic fit requirement [37]. 

(d) Mass window requirements: A mass window re¬ 
quirement on the invariant mass introduces a sys¬ 

tematic uncertainty on its efficiency. The J/ip signal at 
i/s = 4.230 GeV is fitted with a MC shape convoluted 
with a Gaussian function, where the parameters of the 
Gaussian function are left free in the fit. To evaluate the 
systematic effects on the mass window requirement, the 
invariant mass of t^^~ in MC samples are smeared with 
a Gaussian function, where the parameters of the Gaus¬ 
sian function are obtained from the fit. The difference in 
the efficiencies between the signal MC sample with and 







TABLE III. Upper limits of e''"e —>■ tt® J/i/j. The table shows the number of observed events in the tt'^ signal region number 

of events in 7 r° sideband A^.^o, and in J/'i/i sideband efficiency e, the upper limit of signal events with the consideration 

of the selection efficiency fj,~)/e (at the 90% C.L.) and the npper limit of Born cross sections cr^p (at the 90% C.L.). 


^(GeV) 


Ko 

TT^ 


e(%) 

A'^P/e 

(1 + V) 

(1 + A) 

c-^[,(pb) 

3.810 

1 

4 

1 

16.9 

< 20.2 

1.243 

1.056 

< 5.2 

3.900 

0 

1 

2 

29.2 

< 6.0 

0.775 

1.049 

< 2.4 

4.090 

0 

0 

2 

25.7 

< 7.8 

1.078 

1.052 

< 2.2 

4.190 

0 

0 

0 

29.9 

< 6.7 

0.866 

1.056 

< 2.9 

4.210 

1 

1 

1 

29.0 

< 11.8 

0.914 

1.057 

< 3.8 

4.220 

0 

1 

0 

28.5 

< 7.0 

0.937 

1.057 

< 2.2 

4.230 

4 

16 

13 

28.1 

< 18.5 

0.960 

1.056 

< 0.3 

4.245 

1 

1 

2 

27.3 

< 12.6 

0.992 

1.056 

< 3.7 

4.260 

3 

8 

10 

26.5 

< 18.8 

1.021 

1.054 

< 0.4 

4.310 

0 

0 

0 

24.6 

< 8.3 

1.105 

1.053 

< 2.7 

4.360 

2 

3 

4 

23.5 

< 19.9 

1.168 

1.051 

< 0.5 

4.390 

1 

0 

1 

23.1 

< 16.0 

1.198 

1.051 

< 3.9 

4.420 

2 

7 

20 

22.7 

< 16.3 

1.225 

1.053 

< 0.2 

4.470 

0 

0 

4 

22.3 

< 8.9 

1.258 

1.055 

< 1.0 

4.530 

0 

1 

2 

21.8 

< 8.9 

1.295 

1.055 

< 0.9 

4.575 

0 

0 

2 

21.7 

< 9.2 

1.314 

1.055 

< 2.4 

4.600 

3 

5 

7 

21.6 

< 26.2 

1.323 

1.055 

< 0.6 


without mass resolution smearing is 0 . 2 % in the 
mode and 0 . 1 % in the e+e” mode, and is taken as the 
systematic uncertainty. 

(e) Fitting procedure: For the eight data samples with 
clearly observed ry signals, fits to the two photon invari¬ 
ant mass M ( 77 ) are performed to extract the number of 
e+e“ —>■ ryJ/'f/’ decays. The following three aspects are 
considered when evaluating the systematic uncertainty 
associated with the fit procedure. (1) Fitting range: In 
the fit, the M{'y'y) is fitted in a region from 0.2 to 0.9 
GeV/c^. An alternative fit with a different fit range, from 
0.25 to 0.85 GeV/c^, is performed. The differences in the 
yield are treated as the systematic uncertainty from the 
fit range. (2) Signal shape: In the fit, the signal shape 
is described by a shape obtained from a MC simulation 
convoluted with a Gaussian function. An alternative fit 
with a Grystal Ball function [38] for the ry signal shape 
is performed, where the parameters of the Grystal Ball 
function at different CM energies are fixed to those ob¬ 
tained from the fit of the jy signal at ^/s = 4.230 GeV. 
The difference in the yield with respect to the nominal fit 
is considered as the systematic uncertainty from the sig¬ 
nal shape. (3) Background shape: In the fit, background 
shapes are described as a second-order polynomial func¬ 
tion. The fit with a third-order polynomial function for 
the background shape is used to estimate its uncertainty. 
For the data sets where no evident rj signal is found, the 
frequentist method is employed to determine upper limits 
on the Born cross section, and the numbers of signal and 
background events are obtained by counting the entries 
in signal and sideband regions. Two different sideband 
regions, either the ry sideband region or the J/ij: side¬ 
band region, are used to estimate the uncertainty from 
the background shape. The systematic uncertainty as¬ 
sociated with the background shape has been considered 


by taking the most conservative upper limit as the final 
result. 

(f) MUC requirement: In the search for the process 

e“'"e“ —> Tv^J/ip, an additional requirement on the hit 
depth in the MUG for muon tracks was imposed to re¬ 
move the background from e+e“ —>■ By study¬ 
ing the control sample of e“'"e“ —>■ J/ip with a sub¬ 
sequent decay of J/ij) —?► at ^/s = 4.230 GeV, the 

efficiency difference of this requirement between data and 
MG sample was found to be (9.0 ± 1.2)%. The MC effi¬ 
ciency has been corrected for this difference and a value 
of 1 . 2 % is taken as the corresponding systematic uncer¬ 
tainty. 

(g) ISR factor: The uncertainties of the line shape of 
the cross section used in the KKMC generator introduce 
uncertainties in both the radiative correction factor and 
the efficiency. In the nominal results, the line shape of 
the cross section is taken from the fit result from the 
Belle experiment [25]. We have also performed a new fit 
with three incoherent Breit-Wigner functions, including 
the Y (4360) and a second order polynomial function, to 
the same observed cross section a{e^e~ —>■ iqj/ip), where 
the parameters of the Breit-Wigner functions are left free 
in the fit. With this line shape of the cross section, the 
variations in (1 -|- S^) x e are taken as the uncertainties. 

(h) Luminosity: The integrated luminosity of data 
samples used in this analysis are measured using large 
angle Bhabha events, and the corresponding uncertain¬ 
ties are estimated to be 1.0% [39]. 

(i) Branching fractions: The experimental uncertain¬ 
ties in the branching fractions for the processes J/ip ^ 

ry —>■ 77 and 7 r° —>■ 77 are taken from the PDG [34]. 

(j) Other systematic uncertainties: Other sources of 
systematic uncertainties include the trigger efficiency, 
event start time determination and final-state-radiation 
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TABLE IV. Summary of systematic uncertainties (%) in the cross section of e''"e“ 'nJ/'4’ for energies with significant signal 

in the fj,~ {e^ e~) mode. The common uncertainties (Luminosity, Tracking, Photon, Branching fraction and Others) between 
the two modes are shown together. 


Source/v^(GeV) 

4.190 

4.210 

4.220 

4.230 

4.245 

4.260 

4.360 

4.420 

Luminosity 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

Tracking 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

Photon 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

Kinematic fit 

0.4 (0.4) 

0.4 (0.4) 

0.4 (0.3) 

0.4 (0.3) 

0.5 (0.5) 

0.4 (0.4) 

0.3 (0.4) 

0.4 (0.4) 

Mass window 

0.2 (0.1) 

0.2 (0.1) 

0.2 (0.1) 

0.2 (0.1) 

0.2 (0.1) 

0.2 (0.1) 

0.2 (0.1) 

0.2 (0.1) 

Fitting range 

0.0 (1.0) 

0.4 (0.7) 

0.3 (2.6) 

0.1 (2.2) 

0.0 (1.0) 

0.0 (0.6) 

8.6 (7.5) 

0.7 (2.1) 

Signal shape 

0.3 (1.1) 

0.3 (1.1) 

0.3 (1.1) 

0.3 (1.1) 

0.3 (1.1) 

0.3 (1.1) 

0.3 (1.1) 

0.3 (1.1) 

Background shape 

4.6 (0.1) 

3.9 (6.8) 

2.8 (0.0) 

0.0 (0.1) 

9.7 (9.3) 

0.2 (0.0) 

0.5 (0.4) 

0.1 (0.2) 

ISR factor 

0.6 (0.3) 

4.3 (4.2) 

4.7 (6.0) 

4.2 (5.9) 

4.0 (3.6) 

6.6 (5.8) 

9.4 (9.1) 

10.5 (7.7) 

Branching fraction 

0.8 

0.8 

0.8 

0.8 

0.8 

0.8 

0.8 

0.8 

Others 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

Sum 

5.7 (3.6) 

6.7 (8.7) 

6.4 (7.4) 

5.3 (7.2) 

11.0 (10.6) 

7.4 (6.8) 

13.2 (12.3) 

11.0 (8.7) 


TABLE V. Summary of systematic uncertainties (%) in the upper limit on cross section of e'''e iJ/i’ in mode. 


Source/v^(GeV) 

3.810 

3.900 

4.090 

4.310 

4.390 

4.470 

4.530 

4.575 

4.600 

Luminosity 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

Tracking 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

Photon 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

Kinematic fit 

0.4 

0.4 

0.4 

0.4 

0.4 

0.4 

0.4 

0.4 

0.4 

Mass window 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

ISR factor 

0.2 

1.3 

6.3 

0.4 

9.0 

2.8 

1.0 

1.2 

0.8 

Branching fraction 

0.8 

0.8 

0.8 

0.8 

0.8 

0.8 

0.8 

0.8 

0.8 

Others 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

Sum 

3.3 

3.5 

7.1 

3.3 

9.6 

4.3 

3.4 

3.5 

3.4 


simulation. The total systematic uncertainty due to these 
sources is estimated to be less than 1.0%. 

Assuming all of the above systematic uncertainties, 
shown in Table IV, Table V and Table VI, are indepen¬ 
dent, the total systematic uncertainties are obtained by 
adding the individual uncertainties in quadrature. 

For the energy points where statistically significant sig¬ 
nal yields were found, the results from the two J/4’ decay 
modes are found to be consistent. The combined cross 
sections are calculated by considering the correlation of 
uncertainties between these two measurements [40] and 
the results are also listed in Table I. 


VI. SUMMARY AND DISCUSSION 

In summary, using data samples collected with the BE- 
SIII detector at energies from 3.810 to 4.600 GeV, we 
performed an analysis of e+e” —>■ rjJ/tjj. Statistically 
significant r] signals are observed at y/s = 4.190, 4210, 
4220, 4230, 4245, 4260, 4360 and 4420 GeV, and the 
corresponding Born cross sections are measured. In ad¬ 
dition, we searched for the process e+e” —>■ ir^J/ip. No 
significant signals are observed and the upper limits at 
the 90% C.L. on the Born cross section are set. 

A comparison of the Born cross sections a{e'^e~ —>■ 
r]J/il)) in this measurement to that of previous re¬ 


sults [25, 26] is shown in Fig. 4 (a), and a very good 
agreement is achieved. The measured Born cross sec¬ 
tions are also compared to that of e+e” —>■ 'K^'K~J/ip 
obtained from Belle [4] as shown in Fig. 4 (b). Different 
lineshapes are observed in these two processes, which in¬ 
dicate that the production mechanism of the ?7 J/'0 clearly 
differs from that of Tr'^TT~J/ip in the vicinity of ^/s = 
4.1-4.6 GeV. This could indicate the existence of a rich 
spectrum of Y states in this energy region with different 
coupling strengths to the various decay modes. 

The ratio of the Born cross section at 4.260 GeV 
to that at 4.230 GeV, Rt', 23 oi^~^= 
; is calculated to be 0.33±0.04 (common 


systematic uncertainties cancel in the calculation), which 
is found to agree very well with the ratio, —>■ 

wxco) = = 0-43 ± 0.13, of the process 

e+e” —?► wxco [18]. This may indicate that the pro¬ 
duction of riJ/ip and ujXcO are from the same source. 
More data around this energy region may be useful to 
clarify this interpretation. Compared with a theoreti¬ 
cal prediction [24] that considers open charm effects on 
the exclusive cross section line shapes of e+e“ —>■ riJ/tp 
and TT^J/ip, our results on 'qj/ip are within the range of 
the theoretical prediction, and the obtained up¬ 

per limits are higher by a factor of 50 than that of the 
theoretical prediction. More data samples will be helpful 
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TABLE VI. Summary of systematic uncertainties (%) in the cross section of e''"e J/if}. 


Source/v^(GeV) 

3.810 

3.900 

4.090 

4.190 

4.210 

4.220 

4.230 

4.245 

4.260 

4.310 

4.360 

4.390 

4.420 

4.470 

4.530 

4.575 

4.600 

Luminosity 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

Tracking 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

Photon 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

Kinematic fit 

0.4 

0.4 

0.4 

0.4 

0.4 

0.4 

0.4 

0.4 

0.4 

0.4 

0.4 

0.4 

0.4 

0.4 

0.4 

0.4 

0.4 

Mass window 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

MUC cut 

1.2 

1.2 

1.2 

1.2 

1.2 

1.2 

1.2 

1.2 

1.2 

1.2 

1.2 

1.2 

1.2 

1.2 

1.2 

1.2 

1.2 

ISR factor 

0.2 

1.1 

6.5 

0.3 

4.6 

5.7 

3.9 

4.1 

6.7 

0.8 

9.6 

8.7 

7.9 

1.0 

0.7 

0.5 

0.7 

Branching fraction 

0.6 

0.6 

0.6 

0.6 

0.6 

0.6 

0.6 

0.6 

0.6 

0.6 

0.6 

0.6 

0.6 

0.6 

0.6 

0.6 

0.6 

Others 

0.6 

0.6 

0.6 

0.6 

0.6 

0.6 

0.6 

0.6 

0.6 

0.6 

0.6 

0.6 

0.6 

0.6 

0.6 

0.6 

0.6 

Sum 

3.5 

3.6 

7.4 

3.5 

5.8 

6.7 

5.2 

5.4 

7.5 

3.6 

10.2 

9.4 

8.6 

3.6 

3.5 

3.5 

3.5 



\fs (GeV) 



\/s (GeV) 


FIG. 4. A comparison of the measured Born cross section of 
e^e“ —> vJ/i^ to (a) that of a previous measurement [25, 26], 
(b) that of e'’"e“ —>■ TT'^TT~J/'tp from Belle [4]. In these two 
plots, the black square dots are the results of rjJ/'ijj obtained 
in this work and the red star dots are from BESIII(2012). 
The blue dots are results of r)J/ip (a) and tt+tt" J/^/) (b) from 
Belle. The errors are statistical only for Belle’s results, and 
are final combined uncertainties for BESIIFs results. 


to test the predicted cross section of e+e —>■ 
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